The role of single-ion excitations in the mixed-spin 
quasi- l-dimensional quantum antiferromagnet Nd2BaNi05. 

0\ ' A. Zheludev and S. Maslov 






o 
o 



X 



Physics Department, Brookhaven National Laboratory, Upton, New York 11973 USA 

T. Yokoo and J. Akimitsu 
Department of Physics, Aoyama-Gakuin University, 6-16-1, Chitosedai, Setagaya-ku, Tokyo 157 

Japan 



o 

(N 

^', S. Raymond 

Universit de Geneve DPMC 24, quai Ernest Ansermet 1211, Geneve 4, Switzerland. 

o 

,0,: S. E. Nagler 

Oak Ridge National Laboratory, Bid. 7692, MS 6393, P.O. Box 2008, Oak Ridge, TN 37831, 
> 

i^ : USA. 

m 

o 



K. Hirota 
Department of Physics, Tohoku University, Sendai 980-8578, Japan GREST, JST, Japan. 

(February 1, 2008) 

Abstract 

Inelastic neutron scattering experiments on Nd2BaNi05 single-crystals and 
powder samples are used to study the dynamic coupling of l-dimensional 
Haldane-gap excitations in the 8 = 1 Ni^^-chains to local crystal-field tran- 
sitions, associated with the rare earth ions. Substantial interference between 
the two types of excitations is observed even in the l-dimensional param- 
agnetic phase. Despite that, the results provide solid justification for the 



previously proposed "static staggered field" model for i?2BaNi05nickelates. 
The observed behaviour is qualitatively explained by a simple chain-Random- 
Phase-Approximation (chain-RPA) model. 
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I. INTRODUCTION 

Linear-chain rare earth nickelates with the general formula i?2BaNi05 were recently 
recognized as unique model compounds for experimental studies of Haldane-gap quantum 
spin chains in strong staggered applied fields.tl In these materials the Haldane subsystem is 
formed by chains of 5* = 1 Ni^"*" ions with relatively strong, almost isotropic, nearest-neighbor 
antiferromagnetic (AF) exchange interactions (J ~ 27 meV). The corresponding Haldane 
energy gap in the magnetic excitation spectrum is A ^ 11 meV.oTj An effective staggered 
exchange field H^^ acting on these 1-dimensional (1-D) objects is generated by 3-dimensional 
(3-D) magnetic long range ordering in the "auxiliary" sublattice of magnetic rare earth ions, 
that typically occurs at Neel temperatures of several tens of KelvinB The magnitude of the 
staggered field is proportional to the ordered moment on the R^~^ sites, and can be varied 
in an experiment, by changing the temperature. Static magnetic properties of i?2BaNi05 
compounds are rather well described by the chain-Mean Field (MF) model.&Q Using this 
theory to analyze the measured temperature dependencies of magnetic order parameters for 
a number of i?2BaNi05 species, previously enabled us to extract the staggered magnetization 
function for a single Haldane spin chain.tl These experimental results were found to be 
in excellent agreement with theoretical calculations for isolated chains. li^lllil The view of 
i?2BaNi05 systems as Haldane spin chains immersed in a Static Staggered Field (SSF) model 
turned out to be successful in describing the behaviour of Haldane gap excitations as welLM'S 
In agreement with analyticalcJ'Ej and numericaltll predictions, the energy gap was shown to 
increase quadratically with H^^ in Pr2BaNi05, Nd2BaNi05, and (Nda;Yi_3.)2BaNi05 systems 
in their magnetically ordered phases. 

The focus of the present paper is on a different, potentially very important, aspect in the 
magnetism of i?2BaNi05 materials, namely on the strong spin-orbit interactions typically 
associated with rare earth ions. Of particular interest are local (single-ion) crystal field (CF) 
excitations of the rare earths and their possible interactions with the Haldane-gap modes 
propagating on the Ni-chains. Until now this effect has not been investigated in sufficient 



detail. The models used for R^~^ in i?2BaNi05 were in most cases restricted to considering 
the ion's lowest-energy electronic configurations: effective S = 1/2 doublets for Kramers ions 
like Nd^"*" or Er'^+,tl or singlets for Pr'^+ (Ref. 0). Higher-energy electronic states were entirely 
ignored. Yet, in a number of materials, particularly in Nd2BaNi05,Q several CF transitions 
occur very close to the Haldane gap energy. In this situation a strong coupling between 
the two types of excitations can be expected. Does this interaction limit the applicability 
of the SSF model? Does it give rise to peculiar "mixed" Ni-i? excitations? If so, do these 
"mixed" modes retain a predominantly 1-D character of Haldane gap excitations? In the 
present paper we address these important and interesting questions in a new series of inelastic 
neutron scattering measurements on Nd2BaNi05. 

II. EXPERIMENTAL 

The crystal structure of i?2BaNi05 compounds is discussed in detail elsewhere.Ej 
Nd2BaNi05 crystallizes in an orthorhombic unit cell, space group Immm, a ^ 3.8 A, 6 ^ 
5.8 A, and c ~ 11.7 AM The 5* = 1 Ni chains are formed by orthorhombically distorted NiOe 
octahedra that are lined up in the (100) crystallographic direction, sharing their apical oxy- 
gen atoms. The Nd sites are positioned in between these chains, linking nearest-neighbor 
chains along the crystallographic 6-axis (Fig. |l]). We shall denote the corresponding ex- 
change constant as Ji. In addition, one can expect substantial coupling of each Nd'^"*" ion 
to two Ni^+ sites in a third nearby chain (exchange constant J2). The site symmetry for 
the rare earth ions is sufficiently low, and the electronic ground state for Nd^^ is a Kramers 
doublet. The Neel temperature for Nd2BaNi05 is Tn = 48 K. The magnetic structure is as 
described in Ref. p|,p^. 

The simplest technique used to locat CF excitations is inelastic neutron scattering on 
powder samples. Such measurements were previously performed on Nd2BaNi05 for energy 
transfers up to 30 meV using a 3-axis spectrometer.a To learn about higher-energy CF 
transitions we carried out preliminary experiments on the time-of-fiight (TOF) spectrometer 



LAM-D, installed at the spallation neutron source KENS in the High Energy Accelerator 
Research (KEK). LAM-D implements the so-called inversion geometry and is equipped with 
four fixed PG crystal analyzers. The measured spectra are automatically Q-integrated in a 
certain momentum range, that roughly corresponds to 1-2 A~^ at huj = 0. Energy transfers 
from -2 meV to 200 meV are accessible. Energy resolution at the elastic scattering position is 
about 300 /ieV and the relative resolution 5{huj)/Tiuj is fairly constant, around 4%, between 
hu = 2 and 60 meV. In the experiment we used roughly 10 g of Nd2BaNi05 powder prepared 
by the solid state reaction method. The sample environment was a standard closed-cycle 
refrigerator. 

Most of the data presented below were obtained on single- crystal samples. These exper- 
iments were carried out at the High Flux Isotope Reactor (Oak Ridge National Laboratory) 
on the HB-1 triple-axis spectrometer. We used the same assembly of three co-aligned sin- 
gle crystals as in previous studies.!^ The sample was in all cases mounted with the {h, k, 0) 
reciprocal-space plane parallel to the scattering plane of the spectrometer, to compliment 
previous (/i, 0, /) measurements. A neutron beam of 14.7 meV or 30.5 meV fixed final of 
30.5 meV fixed-incident energy was used with 60' — 40' — 40' — 120' coUimations and a py- 
rolitic graphite (PG) filter positioned after the sample. PG (002) reflections were used for 
monochromator and analyzer. The measurements were done with energy transfers up to 
35 meV. The sample environment was a standard Displex refrigerator, and the temperature 
range 35-65 K was explored. The growth directions of the crystals roughly coincide with 
the (100) crystallographic axis. In the experiment the crystal rods were therefore positioned 
horizontally in the scattering plane. In this geometry absorption effects, primarily due to 
the presence of Nd nuclei, may be substantial, as will be discussed in detail below. 

III. RESULTS 



A. Powder data. 

Inelastic powder experiments, while obviously less informative than single- crystal mea- 
surements, provide a convenient view of the excitation spectrum in a wide energy range. 
Figures ^-d show the measured inelastic spectra for Nd2BaNi05 powder at T = 16.5, 55, 
100, 300 K, respectively. The background originating from the sample holder was measured 
separately and subtracted from the data sets shown. Each scan corresponds to 16-20 hours 
of counting time. At T = 55 K, above the Neel temperature Tn = 48 K, one clearly sees 
at least 3 features that are likely to be Nd-CF excitations (Fig. ^). In addition to the two 
previously identified low-energy modes at Ei = 18 meV and E2 = 24 meV energy transfer,Q 
a broad yet intense peak is seen at 40 meV. As in previous 3-axis powder measurements, at 
11 meV one can also see a shoulder corresponding to Haldane excitations. The latter are 
known to have a steep dispersion along the chain- axis. As a result, compared to the dis- 
persionless CF excitations, the Ni-chain modes appear substantially weakened in a powder 
experiment, as opposed to single crystal measurements. No substantial scattering is seen 
within the Haldane gap. At T = 100 K the 11 meV shoulder is much less pronounced, 
while the CF peaks remain almost unchanged (Fig. ^). As could be expected, at high 
temperatures the entire spectrum becomes smeared out (Fig. ^). 

A drastic change in the spectrum occurs upon cooling through the Neel temperature 
(Fig. §a, T = 16.5 K). An intense resolution-limited peak emerges inside the gap and is 
centered at 4 meV energy transfer. This feature corresponds to a dispersion-less Ising-like 
spin wave associated with long-range magnetic ordering.la Below T = 20 K the Haldane gap 
was previously shown to exceed 15 meV.H In apparent contradiction with this single-crystal 
result, the 11 meV shoulder is weakened but still visible in our new 16.5 K powder data. 
One possible explanation of this is that the shoulder is partially due to a weak, previously 
unobserved CF excitation: at low temperatures the Haldane gap increases, but the new CF 
mode stays at 11 meV energy transfer. Some of the single-crystal results presented below 
support this interpretation. 
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A more complete, quantitative, analysis of the powder data is not presently possible, 
as the phonon background is not known. In the future we plan to supplement these data 
with "background" measurements on Y2BaNi05, to isolate the magnetic contribution of the 
rare earth ions in Nd2BaNi05. For the purpose of the present paper however, we can draw 
one useful conclusion: in trying to understand the coupling of the Ni-and Nd-subsystems 
on energy scales comparable to A, we can safely restrict ourselves to considering the two 
strong CF excitations at 18 meV and 24 meV only. 

B. Transverse intensity modulation. 

1. Constant-Q scans 

Having verified the location of major CF excitations in powder measurements, we pro- 
ceeded with more detailed single- crystal experiments. As a first step, we investigated the 
6-axis dependence of intensities and energies for the Ni-chain- and Nd-CF-excitations near 
the 1-D AF zone-center. Typical constant-Q scans collected at Q = (0, A;,2.5), for several 
values of k, are shown in Fig. ^. A total of 9 scans of this type were measured, with k ranging 
from to 1, at T = 55 K (just above the temperature of magnetic ordering, T/v = 48 K). 
These data are summarized in the intensity-contour plot in Fig. |^a. Peaks corresponding 
to the Ni-chain excitations and to the two CF modes are clearly seen. The central result of 
this work is that the relative intensities of these excitations appear to vary significantly with 
transverse momentum transfer. 

To better understand this intensity variation, and to look for a transverse dispersion 
in any of the three modes, we analyzed each of the measured const-Q scans using a semi- 
empirical parameterized model cross section: 

S{Q, oo) = Sh{Q, uj) + SiiQ, oo) + S2{Q, uj), (1) 

A^ 1 

Sh{Q,u;) = Sf^-^-yr^^lf-^J--^2U^M)\ (2) 

S,{Q,u) = Sfh{huj - E,){f^,y{Q))\ (3) 



S2{Q,uj) = Sfh{nu - E2){h^{Q))\ (4) 

{huQf = ^^ + v^sm\Q\\a). (5) 

In these equations Sh{Q,uj), Si{Q,uj), and S2{Q,uj) are the partial contributions of the 
Haldane-gap and two CF modes, respectively. The coefficients Sj^ , Si , and ^2 represent 
the intensities of each component. The spin wave velocity v = 6. OTA |15|and intrinsic energy 
width r characterize the Haldane-gap excitations. In our analysis F was fixed to a value 
of 1 meV, as in previous studies.El Qy is the projection of the scattering vector onto the 
chain axis. Finally, fm{Q) and /Nd(Q) are magnetic form factors for Ni^"*" and Nd^"^ ions. 
The data analysis was performed assuming 7 independent parameters: the three intensity 
prefactors, a constant background, and the three characteristic energies A, Ei and E2. The 
cross section was convoluted with the spectrometer resolution function and fit separately to 
each scan. Typical fits are shown in solid lines in Fig. 0. This analysis yielded a measurement 
of the 6-axis dispersion relation in each branch, plotted with symbols in Fig. ^. To within 
experimental error the CF branches appear flat, and only a very small dispersion is seen in 
the Haldane modes (Fig. ^). The Haldane gap energy is a maximum at the transverse zone- 
boundary k = 0.5. The intensity variation for each mode is shown in Fig. ^. The intensity 
of the 18 meV excitation appears almost /c-independent. At the same time, intensities of 
both the Haldane- and the 24 meV CF- modes are strongly fc-dependent, and are seemingly 
in "antiphase" one to another. 

Data sets similar to those described above were also collected at T = 35 K, i.e., below 
the Neel temperature. These results are summarized in Figs. ^ and |^. Compared to the 
high-temperature measurements, the Haldane gap energy is increased to about 14 meV. The 
Haldane branch thus occurs closer to the 18 meV CF mode. It is also almost twice as weak as 
at T = 55 K. As a result, the scattering of data points in the measured transverse dispersion 
curve shown in Fig. ^ is rather large, and the bandwidth can not be reliably determined. 
Compared to the high-temperature regime, the magnitude of intensity variations in the both 
the Haldane- and 24 meV- excitations are visibly reduced (see Fig. |^). 



2. Global fits to data 

The final analysis of the measured constant-Q scans was performed by fitting a model 
cross section simultaneously to all data collected at each temperature (Figs, ^a and ^). 
The formulas [^-|^ used to analyze individual scans were modified to allow for a transverse 
dispersion and intensity modulation in each branch. The appropriate forms for these k- 
dependencies are derived in the theory section below. The corresponding expression for the 
magnetic dynamic structure factor near the 1-D AF zone-center is as follows: 

S{Q,uj) = Sh{Q,uj) + S^{Q,uj) + S^{Q,uj), (6) 



:(o) 



A^ 1 



r {hujy + A'^{huj - huQY/v^) 



Sh{Q,u;) = S';;' — -—^——- , ,,,^,, {fm{Q)r + AH{hi{Q)fm{Q))cos{Q^b/2) , (7) 



.(0); 



Si{Q,uj) = s\'"6{nuj~E,{Q^)) (/Nd(Q)r + Ai(/Ni(g)/Nd(g))cos(gx&/2) , (s) 



.(0); 



S2{Q, uj) = S^S {huj - E^iQ^)) {hd{Q)y + ^2(/Ni(g)/Nd(g)) cos{Q^b/2) , (9) 



(hiUQ) = ^A2 + t;2 sin2(g||a) + Bh cos^iQ±b/2), (10) 

Ei{Q^) = Ei + Bicos\Q^b/2), (11) 

E2{Q±) = E2 + B2 cos\Q^b/2). (12) 

In these equations Q± is the component of momentum transfer parallel to the b axis. As 
explained in the theory section, the dimensionless coefficients Ah, A\ and A2 quantify the 
admixture of Nd-CF spin fluctuations to the chain modes and vice versa. Bh, B^ and B2 
are respective magnitudes of transverse dispersion. To minimize the number of adjustable 
parameters the values -Bi = 0, ^2 = (no dispersion observed) ,and Ai = (the mode 
appears not to be modulated) were fixed, which hardly affected the final x^- The results of 
the global least-squares fitting to the T = 55 K and T = 35 K data are shown in Figs. ^3 
and ^, respectively. The refined values of all parameters are summarized in Table |. In 
both cases a x^ = 1-7 was achieved, which suggests that the model cross section describes 
the observed behavior reasonably well. As discussed above, the coefficients for intensity 
modulations in all three modes, as well as the dispersion bandwidth for the Haldane branch 
at T = 55 K, determined in this analysis, can be considered rather reliable. On the contrary. 



the large value for Sh at T = 35 K is somewhat suspect. In the refinement this parameter 
has a strong correlation with Ei, i.e., the Haldane gap is not very well resolved. This also 
accounts for the large scattering in the points obtained in fits to individual scans in Fig. |^. 
There is little doubt however that the bandwidth of transverse dispersion remains very small 
even in the magnetically ordered state. 

3. Wide-range constant-E scans 

It is implied in Eqs. ^-^ that the intensity modulation in all modes is periodic along the 
6-axis. To verify this periodicity we performed constant-ii^ scans at the Haldane gap energy 
and at huj = 25 meV, covering a broad range of transverse momentum transfers. Typical raw 
data collected at T = 55 K are shown in Figs. ^ and p!0|a. A serious technical problem with 
such measurements is that they are strongly infiuenced by neutron absorbtion in the sample. 
In each scan k varies significantly at a constant value of h. As a result, the crystal rods 
rotate relative to the incident and outdoing beams by a large angle. Absorbtion corrections 
could be calculated, but additional undesirable effects arising from the size (length) of the 
sample being comparable to the width of the incident neutron beam, are difficult to account 
for analytically. To correct the data for any transmission effects, we measured the actual 
transmission coefficient in situ. This was done separately for each point of the constant-i? 
scans performed. The monochromator, sample and scattering angles of the spectrometer 
were set to values calculated for a given point in the inelastic scan. The analyzer angles 
were then set to the elastic position to measure the intensity of incoherent scattering. Inco- 
herent scattering being isotropic, the measured incoherent intensity can be assumed to be 
proportional to the effective transmission coefficient. Typical measured transmission curves 
are shown in Figs. ^ and p^. A few Bragg reflections were accidentally picked up in these 
elastic scans, and the corresponding points were removed from the data sets shown. In 
Figs. ^3 and p!0[ b arrows indicate geometries of minimum transmission, realized when the 
sample rods are parallel to the incident or the scattered beams. The measured transmission 
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coefficients were used to scale the measured const-i? scans point-by-point. The resulting 



corrected data are plotted in Figs. and 11 



Figure ^ shows transmission-corrected /c-scans taken at the Haldane gap energy, near two 
1-D AF zone-centers h = 1.5 and h = 2.5. The intensity fall-off at large \k\ is mostly due to 
focusing and to the effect of magnetic form-factors. On top of this fall-off, a periodic intensity 
modulation is clearly seen. The data can be rather well fit to Eq. |^ convoluted with the 
spectrometer resolution function. The fits are shown in solid lines in Fig. ^ and correspond 
to An = —0.19(0.02). This value is in good agreement with A-^ = —0.17, obtained from the 
analysis of constant-Q scans described above. 



Figure pT| b shows that near the 1-D AF zone-center the intensity of the 25 meV CF 
mode is modulated as well, in antiphase with that of the Haldane branch. The data can be 
fit to Eq. ^ assuming A2 = 0.28(0.05) (solid line in Fig. |ll|b). Note that this value is less 
than a half of that deduced from the analysis of constant-Q scans at the same temperature 
{A2 = 0.72). This discrepancy is easily understood. Even at 25 meV energy transfer the 
Haldane-gap excitations contribute to inelastic scattering, thanks to their steep dispersion 
along the chain axis that produces long "tails" on the high-energy side of the peak in 
constant-Q scans. The contribution of these tails was not taken into account in the analysis 



of the wide-range const-E' scan shown in Figure [Tl|b, but was included in the global fit to 
constant-Q scans. Our previous estimate for A2 should thus be considered more accurate. 

An important result is that away from the 1-D AF zone center, at h = 2.25, the 25 meV 
mode shows no oscillations of intensity (Fig. |ll]a). To within experimental error the measured 
k dependence can be explained by the Nd-form factor alone (solid line). This observation 
supports the notion that intensity oscillations result from a mixing between CF and Haldane- 
gap excitations. Away from the 1-D AF zone-center Ni-chain modes occur at very high 
energy transfers and therefore do not interact with the dispersionless CF branch. 
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4- Previous data analysis procedures revised 

The discovery of a substantial variation of the intensity in the CF excitations raises 
concerns about the procedure used to analyze the inelastic data for Nd2BaNi05 in Refs. ^,|^. 
In these studies the CF "background" was measured away from the 1-D AF zone-centers, 
where the Haldane modes are not seen. The background was then subtracted point-by- 
point from the "signal" scans measured at the 1-D AF zone-center. It was assumed that 
such background-subtracted scans are a good measure of inelastic scattering originating 
from the Ni-chains alone. We now understand that this approach is, in principle, not valid. 
Although CF excitations are almost dispersionless, their intensity in the "background" scans 
will not necessarily match that in the "signal" scans. According to the theoretical model 
described below, the procedure can only be applied at A; = {2n + l)/2, h = (2m + l)/2 (m, 
n- integer), where there is no mixing between Ni-chain and Nd-single-ion excitations. The 
data analysis employed in Ref. ^, where measurements were performed at Q = (1.5,0.5,0), 
is thus fully appropriate. In Ref. ||, on the other hand, the data were taken at k = 0, and 
the background subtraction may, in principle, have produced systematic errors. It appears 
however that point-by-point background subtraction is a fairly reliable way to determine 
the Haldane gap energy at any transverse wave vector. Indeed, in Ref. ^ all scans were 
taken below 25 meV energy transfer, and thus include only the 18 meV CF excitation. This 
mode, as we now know, is in fact fiat. The main results of Ref. ^ pertain to the temperature 
dependence of the Haldane gap, that occurs even below 18 meV in the studies temperature 
range, are totally valid. 

C. Dispersion along the chain axis 

The mixing between the Haldane-gap and 24 meV CF modes should be most prominent 
at points of reciprocal space where the two excitations have similar energies. While the 
CF excitation is practically dispersionless, the Haldane branch has a steep dispersion along 
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the chain axis, and can be expected to anticross with the 24 meV branch rather close 
to the 1-D AF zone-center. To investigate this phenomenon we measured the intensity 
of inelastic neutron scattering from Nd2BaNi05 on a grid of points m h — E space at 
k = 0.5, h = 1.5 — 1.8, at T = 55 K. The result of this measurement is summarized in 
the grayscale/contour plot in Fig. O. The circles and squares show the positions of peaks 



seen when the data are broken up into a collection of constant-Q and constant-i? scans, 
respectively. The most prominent feature of the spectrum is an anticrossing at the point of 
intersection of the flat 24 meV CF mode and the parabolic Haldane branch. The solid lines 
are guides for the eye that emphasize this effect. The anticrossing gap 26 E can be estimated 
to be roughly 4 meV. Note that there is no evidence of anticrossing of the Haldane and 
18 meV-CF branches. 

A much less significant feature that is barely seen within the statistics of the measure- 
ment is what appears to be a new dispersionless mode at hu ^ 12 meV (dash-dot line in 
Fig. |T2|). In the single crystal sample it appears much weaker than the the Haldane mode, 
but nevertheless may account for the "shoulder" seen in powder samples at T = 16.75 K, 
as discussed previously. This new feature is substantially weaker than the other two CF 
excitations, and we shall disregard it in the following analysis. 

IV. THEORY 

It is rather interesting that in Nd2BaNi05 the pronounced transverse modulation of in- 
tensity in the Ni-chain and CF excitations is not accompanied by any significant transverse 
dispersion. If we were dealing with a one-component system, e.g., directly coupled chains, a 
narrow bandwidth along a given direction would automatically result in only a weak period- 
icity in the structure factor. For magnetic excitations in a Bravais crystal, this follows from 
the sum rule for the first moment of the magnetic dynamic structure factor S{q,uj)J^ The 
strong modulation and lack of dispersion in Nd2BaNi05 must therefore be an interference 
effect, that involves both Ni- and Nd magnetic degrees of freedom. In the following sec- 
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tions we shall demonstrate that this mixing can be qualitatively accounted for by a simple 
Random Phase Approximation (RPA) model for Ni-R interactions. 

A. Bare susceptibilities 

1. Ni- chains 

The first step in the RPA analysis is to determine the bare (non-interacting) dynamic 
susceptibilities for the Ni and Nd sublattices. For the Haldane spin chains in the vicinity of 
the 1-D AF zone-center q\\ = vr we can safely use the single mode approximationcJll^: 

(o)/r> ^ _ l-cosQ||a Zv 

In this formula Q\\ is the projection of the momentum transfer Q onto the chain axis, v is 
the spin wave velocity in the chains, v ~ 6. OTA, and Z ^ 1.26 (Ref. [T5| ). It is easily verified 
that for scattering vectors close to the 1-D AF zone-center Eq. |^ used to fit the the inelastic 



scans coincides with the imaginary part of Eq. 13 in the limit e ^ 0. In the magnetically 
ordered state, i. e., in the presence of an effective mean field Heff, Eq. (|13[) can still be 
utilized. However, the corresponding gap energies for the three components of the Haldane 
triplet will be modified by H^fj, as discussed in Refs. P|, pi]| , |2D| J5|. 



2. Rare earth ions 

For the rare earth subsystem the bare single-ion susceptibility can not be written down 
so easily, as the exact electronic configuration of i?^^ in i?2BaNi05 has not been calculated 
to date. When we were concerned with the mechanism of magnetic ordering in Nd2BaNi05 
in Refs. HH, we were dealing with the static susceptibility of the rare earths. For Nd'^"^ at 
low temperatures the static susceptibility can be, to a good approximation, attributed to 
the lowest-energy Kramers doublet and written as a S = 1/2 Brillouin function. For our 
present purpose of studying the dynamic coupling of Nd and Ni subsystems we are interested 
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in the dynamic susceptibility of Nd^"*" at frequencies comparable to the Haldane gap energy. 
In the paramagnetic phase the contribution to XNd from the ground state doublet is purely 
elastic and can thus be totally ignored in the RPA calculation. At T > Tn we therefore only 
need to consider the contribution of transitions between the ground state of the ion and 
higher-energy crystal-field levels. In the ordered phase, we also have to take into account 
single-ion transitions within the ground state double, which becomes split by the mean 
exchange fields. As a first step, we shall consider only one excited level of Nd^^, at energy 
E. Transitions from the ground state to the excited state produce the following term in Nd 
single-ion susceptibility: 

X^Zi^) = M^ ^._^l^^^y_ R(T) (14) 

In this formula M is the matrix element of the magnetic moment operator between the 
ground state and the excited state. R{T) a the temperature renormalization factor defined 
as the difference in populations of the excited and ground states. To simplify our qualitative 
analysis, in Eq. |T^ we have ignored any terms that are off-diagonal in spin projection indexes, 
and the above formula thus applies to a particular selected channel of spin polarization. 

B. Interaction geometry 

Let us now consider the geometry of Ni-Nd exchange interactions (Fig. [l|). The an- 
tiferromagnetic coupling denoted by Ji stabilizes the magnetic structure seen in the low- 
temperature phase, and is responsible for the staggered mean field at T < Tn. As is often 
the case in body-centered crystals, the coupling J2 is geometrically frustrated. Indeed, it 
links a single R^^ moment to two subsequent Ni^"*" spins in the chains, that are strongly 
coupled antiferromagnetically between themselves. As a result, at the MF level J2 plays 
no role in long-range ordering. Below we shall arbitrarily assume that the exchange matrix 
is diagonal in spin projection indexes. In the body-centered structure, with 2 i?-sites per 
every Ni, we have to consider 6 magnetic ions: Ni(l) at (0.5,0.5,0), Ni(2) at (0,0,0.5), Nd(3) 
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at (0.5,0,d), Nd(4) at (0.5,0,-d), Nd(5) at (0,0.5,d-0.5), and Nd(6) at (0,0.5,-[d-0.5]), where 
d ~ 0.2. For the RPA calculation we need the Fourier transform of the exchange matrix in 
this basis set: 



/ 



JdQ) 



Ji(Q) Ji{Q) J2(Q) jm)^ 

J2(Q) r^{Q) JiiQ) J*AQ) 

JiiQ) ^2*(Q) 
MQ) J2{Q) 

j;{Q) JiiQ) 

J2{Q) MQ) ^ 

Ji{Q) = 2Ji cos(7rA;) exp(27ri/(i) 
J2{Q) = 2J2 cos(7r/;,) exp[27Til{d - 0.5)] 



(15) 



(16) 
(17) 



One readily sees that, in the reciprocal-space planes {h, ^ ^ (^-integer), Ji-interactions 



cancel out and can be ignored. A similar cancellation occurs for J2 in the 



'2n+l 



,k,l) (n- 



integer) planes. These cancellations result from interaction topology, and do not rely on any 
of our assumptions regarding the bare susceptibilities. 

C. RPA susceptibility 

The assumptions made above bring us to a simple mode-coupling problem in each 
channel of spin polarization. The RPA susceptibility matrix is written as x^^^(Q)^) = 
[l + x'{Q,uj)J{Q)r'x%Q,u;). 



1. Dispersion 



RPA 



0. In our 



The dispersion relation is obtained by solving the equation det I x 
case of only one excited Nd-state this can easily be done analytically. With 6 independent 
atoms we obtain 6 modes (two of them degenerate): 
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2{hui,2? = A^ +E^± ./(A2 - ^2)2 + sEM^ZvJl 



hiu^fi = E, (18) 

where we have introduced the notations Ag = A^ + f ^ sin^(Q||a), and 
J| = I Ji(Q) ±J2{Q)f = AJl cos'^ink) + 4J| cos^ (Trh) ± 16Ji J2 cos(7rA;) cos(7r/i) cos(7r/). 

(19) 
We have imphcitly dropped the (1 — cosQ||a)/2 factor in Xm ^^d the temperature factor R 

(0) 

m XNd- 

Particularly useful is the linearized version of Eq. |l^ which becomes valid in the weak 
coupling limit: 

, Zv EM^ 



hu2{Cl) = E + J+{Qf 



Q\\ 
Zv EM^ 



2 ' 



A ^^-A^^^ 



A E^-A% ' 



huj,{Cl) = E + J4Qf^'' ^^' 



^E^- ^Q, ' 



hoo^AQ) = E, (20) 

The 1st and 3rd branches can be described as Ni-chain Haldane gap modes with an admixture 
of Nd-spin correlations. The 2nd and 4th branches, conversely, are CF excitations with a 
small component of Ni-spins. If the Nd-single-ion excitation lies above the Haldane branch at 
a given wave vector, the energy of the Ni-chains mode is pushed down by Ni-Nd interactions, 
while the CF excitation energy is increased. The reverse is expected in the case of a single-ion 
excitation inside the Haldane gap. Modes 5 and 6 appear as totally unperturbed single-ion 
transitions. A very important result is that the transverse dispersion resulting from Ni-i? 
interactions is quadratic in J' and is thus expected to be very small in the weak-coupling 
limit. 
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Also useful for our purposes are values of anticrossing gaps at the point of mode inter- 
section. According to Eqs. [13 these are given by: 



25Eu 



J+M\ 



l2Zv 



25^3^4 = J-mJ'^, 



(21) 



for the two pairs of intersecting same-symmetry branches, respectively. Intersection of 
branches of different symmetry (1,4) and (2,3) does not lead to an anticrossing effect. 



2. Intensity modulation 

The inelastic neutron scattering cross section is related to the imaginary part of dynamic 
susceptibility through the fluctuation-dissipation theorem. To get the actual intensities we 
also have to take into account the Ni- and R- magnetic form factors /Ni(<5) and f-^diQ)'- 

^ hi{Q) ^ 



da 



dn 



■IF^'^ fm{Q) fmiQ) fm{Q) fmiQ) fmiQ) /Nd(Q) )x"{Q,^ 



V 



22) 



/ 



/Ni(Q) 

/Nd(g) 

fmiQ) 

fmiQ) 
fmiQ) 

The task of calculating the total x"iQ,^) is greatly simplified in the weak-coupling limit, 
where one can write x^^^iQ^^) ~ [1 ~ X^iQj^)JiQ)]x^iQj^)- At this level of approx- 
imation the two Haldane excitations (modes 1 and 3) are not distinguished. The same 
applies to CF modes 2,4,5,6. For the sum of partial cross section of the Haldane modes a 
straightforward yet somewhat tedious calculation gives: 



da 



1+3 



Zv 



dfl 



oc 



A, 



/Ni(g)'-4/Ni(g)/Nd(Q) 



EM^ 



E^-A% 



-Re[MQ) + J2iQ)]\. 



(23) 



Similarly, for the single-ion modes: 



^ oc 2M2 [ /Nd(Q)' + 4/Ni(Q)/Nd(Q)-^^^^Re IMQ) + UQ)] I . (24) 

One readily sees that, unlike the corrections to the dispersion relation, corrections to mode 
intensities are linear in Ji and J2 and can be substantial. The linear form of Eqs. pD|,pB| and 



23 enables a straightforward generalization to the case of several Nd-single-ion transitions. 



3. Results for the 1-D AF zone-center 

The expressions for the dynamic structure factor can be further simplified for scattering 
vectors in ( ^"^^ , fc, /) (n-integer) reciprocal-space planes, where most of the neutron data 
were collected. Here the effect of J2 is totally cancelled, and there are only 3 branches in the 
spectrum: the Haldane branch (labeled "H"), a CF branch that interacts with the Ni-chain 
excitations ("CF"), and the unperturbed CF transition ("CFq"). Using Eq. ^on the exact 
(rather than linearized) RPA susceptibility matrix, and truncating the result to the first 
order in Ji we obtain: 



^ oc 1^ ih,{Qf - 8J^hi{Q)fMQ)^^, 



« ^ I hiiQ? - 8Jihi{Q)fMQ) ^,^^'\, cos(7rfc) cos(27rM) ) , (25) 

Oil \ ^ -^Qii 



^oc2M2 [/Nd(Q)2[l + cos(47r/rf)] + 8Ji/Ni(Q)/Nd(Q)^^^^^cos(7rA;)cos(27r/rf)J , (26) 

^ oc 2M2/Nd(Q)'[l - cos(47r/rf)]. (27) 

The intensity of all 3 modes become modulated. The period of modulation along the 6-axis is 
exactly 26*, while that along c is incommensurate with the lattice. For the reciprocal-space 
rods ( ^"^^ , k, 0) (n-integer) the mode CFq vanishes and only 2 branches remain. For the 
purpose of convenience we shall also re- write Eqs. ^ for the (^^^, k, I) (n-integer) planes: 

;.^h(Q) = Aq, - 4Ji(Q)^^ ^f_^^, cos'ink), (28) 

^^cf(Q) = Ai+4Ji(Q)^^ ^f_^^, cos^(7rfc), (29) 

^^CFo(Q)=^- (30) 
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As noted previously, the magnitude of transverse dispersion is proportional to Jf. 

We can not expect the model described above to be quantitatively valid for any real 
i?2BaNi05 system. Its main limitation is that it assumes the exchange matrix and the 
matrix elements of momentum between the rare earth ground- and excited states to be 
simultaneously diagonal in spin projection indexes. In the i?2BaNi05 structure, the rare 
earth sites have very low symmetry and the bare susceptibility and exchange matrixes are 
bound to be rather complex. We can however expect the model to make qualitatively 
correct predictions regarding the periods and signs of modulations of excitation intensities 
and energies. 

4- Paramagnetic phase 

For Nd2BaNi05 we should distinguish two regimes. At T > Tn, at energy transfers 
below 30 meV, only the 18 meV and the 24 meV Nd-modes will mix with the Haldane-gap 
excitations. Experimentally, the 18 meV transition appears totally decoupled from the Ni- 
chains . This is most likely due to its particular polarization and the actual form of the 
exchange tensors. For example, if the 18 meV excitation were polarized along the x axis, 
and Jf^ is zero, than no mixing would occur. The 24 meV excitation is thus the only one 
that mixes with the Haldane branches, producing the observed modulations. 

The RPA model, oversimplified as it is, justifies the use of Eqs. |6|-[I^ in fitting the inelastic 



neutron scattering data. From this analysis we can even extract some useful quantitative 
information. According to Eqs. ^ and ^ the ratio Bn/An (see Eqs. [7|,p!0D is independent 
of any parameters of the rare earth ions: 

Zv 
Bn/An = — Ji (31) 

This is a robust result that allows us to obtain at least an order-of-magnitude estimate for 
Ji, without making any assumptions regarding Nd^"*". For an isolated Haldane spin chain 
Zv ~ 7.6A, so Bn/An ^ 3.8Ji. From the T = 55 K experimental values A = —0.17 meV 
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and B = —0.19 we get Ji ~ 0.3 meV. We can compare this value to the MF couphng a 
for Nd2BaNi05 that we obtained in Refs. y,^ Indeed, Jj_ should be related to a through 
JiM^^^^ M^^'^'' = a/2, which gives Ji ~ 0.5 meV, in reasonable agreement with our new 
estimate. 

A more crude estimate can also be obtained for J2, which in our model is the only 
possible source of the observed anticrossing between the Haldane and 24 meV CF branches 
(Fig. |12]). Indeed, at A; = 0.5 where these data were taken the effect of Ji is cancelled, 
according to Eq. [T^. Experimentally, the anticrossing is about 4 meV. From Eq. ^ we 
get J2M1 ^ 0.8 meV. The 24 meV CF mode and the Haldane excitation are of comparable 
intensity, and based on Eqs. ^ and ^ we can assume that M^ ~ Zv/A, which gives M ~ 3. 



This gives us an estimate for J2: of the order of 0.2 meV. This is a rather reasonable value 
for a rare earth- transition metal superexchange bond. 

5. Ordered state 

In the magnetically ordered phase an additional Nd-centered excitation, namely the 
4 meV Ising-like spin wave that appears below T^ (Fig. ^ja) will also mix with the Ni-chain 
modes. Being inside the Haldane gap, it will produce intensity and energy modulations 
that are opposite to those resulting from the mixing with the high-energy 24 meV branch. 
Intensity oscillations in the Haldane branch will thus be reduced, in agreement with our ex- 
perimental findings. As observed in previous studies, the intensity of the peak corresponding 
to Haldane excitations is reduced by almost a factor of 2. This results from the suppression 
of the 6-axis-polarized Haldane excitation (longitudinal mode) in the ordered state.Q In the 
framework of the RPA model, a decrease of spectral weight in the Haldane branches in- 
evitably leads to a smaller interference effect in the CF mode. Such a reduction of intensity 
modulation in the ordered state is indeed observed experimentally for the 24 meV excitation 
(see Table %. 
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V. SUMMARY 

We have demonstrated that local excitations associated with magnetic rare earth ions 
play an important role in the magnetism of i?2BaNi05 quantum antiferromagnets. Even in 
the 1-D (paramagnetic) phase there are substantial dynamic interactions between quantum 
spin chains and the rare earth subsystem. This phenomenon can be qualitatively accounted 
for by a simple chain-RPA model. Not having more detailed information on the electronic 
states of the rare earth ions, using this model for a quantitative analysis of the data is 
obviously a leap of faith. It does however produce reasonable order-of-magnitude estimates 
for Ni-Nd coupling constants. 

At any temperature the Haldane-gap modes can not be considered as purely-Ni exci- 
tations, and their interference with Nd-single-ion modes is substantial. Despite that, the 
Haldane gap energy is practically the same as in isolated chains in a static mean exchange 
field, and the transverse dispersion is very small. The Static Staggered Exchange Field model 
is thus fully justified for Nd2BaNi05, as are all our previous results for the temperature- 
(staggered field-) dependence of the Haldane gap. 
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FIGURES 
FIG. 1. A schematic representation of the magnetic interaction geometry in the structure of 

i?2BaNi05. 

FIG. 2. LAM-D inelastic spectra measured in Nd2BaNi05 powder samples at different temper- 
atures. The contribution of the sample holder has been subtracted. 

FIG. 3. Typical constant-Q scans measured in Nd2BaNi05 at T = 55 K at the 1-D AF 
zone-center h = 2.5, at different momentum transfers perpendicular to the chain axis (symbols). 
The solid lines are a global fit to the model cross section, as described in the text. The hatched, 
light-gray and gray peaks represent partial contributions from the Haldane-gap and two crystal-field 
excitations, respectively. 

FIG. 4. (a) Contour plot of inelastic scattering intensity measured in Nd2BaNi05 at T = 55 K 
at the 1-D AF zone-center h = 2.5 in a series of constant-Q scans (20 counts/3 min. contour steps). 
The symbols show the energies of Haldane (circles), and two crystal- field excitations (triangles), 
as determined from fitting a model cross section to individual scans, (b) The result of a global fit 
to the data shown in (a). 

FIG. 5. (a) A;-dependence of excitation intensities measured in Nd2BaNiOl5 at T = 55 K, as 
deduced from model cross section fits to individual const-Q scans (see Fig. |^). (b) fc-dependence 
of Haldane gap energy, determined in the same fashion. Solid lines are guides for the eye. 

FIG. 6. (a) Contour plot of inelastic scattering intensity measured in Nd2BaNi05 at T = 35 K 
at the 1-D AF zone-center h = 1.5 in a series of constant-Q scans (20 counts/3 min. contour steps). 
The symbols show the energies of Haldane (circles), and two crystal- field excitations (triangles), 
as determined from fitting a model cross section to individual scans, (b) The result of a global fit 
to the data shown in (a). 
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FIG. 7. (a) ^-dependence of excitation intensities measured in Nd2BaNi05 at T = 35 K, as 
deduced from model cross section fits to individual const-Q scans (see Fig. |^). (b) fc-dependence 
of Haldane gap energy, determined in the same fashion. Solid lines are guides for the eye. 

FIG. 8. Measuring the transverse modulation of the Haldane-gap excitations intensity: (a) 
Constant-i? scan measured in Nd2BaNi05 at T = 55 K near the 1-D AF zone-center h = 1.5 at 
huj = 11.5 meV (open circles). The background is measured off the 1-D AF zone-center (solid 
circles), (b) Transimission curve for the scan shown in panel (a), measured as explained in the 
text. 

FIG. 9. Transverse modulation of the Haldane-gap excitations intensity: Constant- i? scans 
measured in Nd2BaNi05 at T = 55 K near the 1-D AF zone-centers h = 1.5 (a) and h = 2.5 
(b) at huj = 11.5 meV (symbnols). The measured intensity has been corrected for absorbtion and 
the measured background has been subtracted. Solid lines are fits to a model cross section, as 
described in the text. 

FIG. 10. Measuring the transverse modulation of the 25 meV excitation intensity: (a) 
Constant-£' scan measured in Nd2BaNi05 at T = 55 K near the 1-D AF zone-center h = 2.5 
at huj = 25 meV. (b) Transmission curve for the scan shown in panel (a), measured as explained 
in the text. 

FIG. 11. Transverse modulation of the 25 meV excitation intensity: Constant- -E scans mea- 
sured in Nd2BaNi05 at T = 55 K away from (a) and close to (b) the 1-D AF zone-center h = 2.5 
at huj = 25 meV (symbols). The measured intensity has been corrected for absorbtion and the 
measured background has been subtracted. Solid lines are fits to a model cross section, as described 
in the text. 
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FIG. 12. Logarithmic contour plot of inelastic neutron scattering intensity measured in 
Nd2BaNi05 at T = 55 K near the 1-D AF zone-center h = 1.5. Circles and squares show peak 
positions as they appear in constant-i? and constant-Q scans, respectively. The bars are peak 
widths at half height. Solids lines are guides for the eye. Dazhed lines represent the dispersion 
relation expected for the Haldane excitations (parabola) and 25 meV CF mode (straight line) in 
the absense of mode coupling. The dash-dot line shows the position of what seems to be a weak, 
previously unobserved CF excitation. 
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TABLES 
TABLE L Parameters characterizing the transverse dispersion and intensity modulation of the 

Haldane-gap and CrystaLfield excitations in Nd2BaNi05, as obtained from global fits of a model 

cross section to series of constant-Q scans. 

T = 55 K Q = (2.5, A;, 0) T = 35 K Q = (1.52, k, 0) 

A (meV) 10.97 (0.09) 13.26 (0.07) 

An (meV) -0.17 (0.01) -0.12 (0.02) 

Sh -0.19 (0.14) 0.64 (0.15) 

S^^ (arb. u.) 1.17 (0.03) 0.66 (0.02) 

El (meV) 18.33 (0.09) 19.7 (0.1) 

Ai (meV) (fixed) (fixed) 

Bi (fixed) (fixed) 

S[^^ (arb. u.) 0.20 (0.01) 0.22 (0.02) 

E2 (meV) 23.8 (0.06) 24.6 (0.09) 

A2 (meV) (fixed) (fixed) 

B2 0.72 (0.06) 0.37 (0.05) 

4°^ (arb. u.) 0.34 (0.01) 0.30 (0.01) 
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